New classical brackets for dissipative systems. 
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A set of brackets for classical dissipative systems, subject to external random forces, are derived. 
The method is inspired to the old procedure found by Peierls, for deriving the canonical brackets 
of conservative systems, starting from an action principle. It is found that an adaptation of Peierls' 
method is applicable also to dissipative systems, when the friction term can be described by a linear 
functional of the coordinates, as is the case in the classical Langevin equation, with an arbitrary 
memory function. The general expression for the brackets satisfied by the coordinates, as well as 
by the external random forces, at different times, is determined, and it turns out that they all 
satisfy the Jacobi identity. Upon quantization, these classical brackets are found to coincide with 
the commutation rules for the quantum Langevin equation, that have been obtained in the past, by 
appealing to microscopic conservative quantum models for the friction mechanism. 
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I. INTRODUCTION 

The study of quantum dissipative systems, has at- 
tracted, in the last decades, a lot of interest, in view of 
its broad spectrum of applications, ranging from quan- 
tum optics through statistical mechanics, etc. The stan- 
dard approach to deal with quantum dissipation, is based 
on the idea that the physical origin of dissipation is the 
interaction of the system with a heat bath, consisting 
of a large number of degrees of freedom. One considers 
then some microscopic, conservative model for the heat 
bath (and its interaction with the system), and tries to 
recover the macroscopic quantum behavior of the dissipa- 
tive system alone, by eliminating from the description the 
degrees of freedom describing the bath. In Ref.(|lJ), it is 
shown, indeed, that the most general quantum Langevin 
Equation, which is one of the most popular models for 
dissipation, can be obtained from a simple microscopic 
model, where the heat bath is described by a set of in- 
dependent oscillators, linearly coupled to the system of 
interest. 

One may wonder whether it is possible to find a quan- 
tization method for dissipative systems, which is based 
on the macroscopic description of dissipation only, and 
makes therefore no use of microscopic models. As is 
well known, quantization of dissipative systems is by 
no means straightforward, because in general they ad- 
mit neither a Lagrangian nor a Hamiltonian formulation. 
Moreover, even in those special instances where such a 
formulation can be given, the application of the conven- 
tional canonical quantization methods leads to physically 
incorrect results 0- In this Letter, we show that new 
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classical brackets can be consistently built for dissipative 
systems, by generalizing the covariant definition of Pois- 
son Brackets for Lagrangian systems, discovered long ago 
by Peierls (See also Refs.0, H H 0). Our bracket is 
defined on the infinite-dimensional functional space con- 
sisting of all possible classical trajectories, that are ac- 
cessible to the system under the influence of the random 
force. It turns out that, when dissipation is present, the 
random external force has a non-vanishing bracket with 
the system coordinates, which implies that it cannot be 
consistently taken to be zero. This seems to be in agree- 
ment with the fluctuation-dissipation theorem, which re- 
quires fluctuating forces, in the presence of dissipation. 

By the correspondence principle, our classical brackets 
can be eventually quantized, upon substituting them by 
(l/(ih) times) commutators. In this way, we recover the 
same expressions for the commutators between the sys- 
tem coordinates and the random forces, which were de- 
rived from the independent oscillator microscopic model 
of Ref.Q. 

In what follows, we make no attempt at mathematical 
rigor, and the presentation is totally heuristic. We hope 
to clarify elsewhere the delicate issues involved in the con- 
sideration of Poisson structures in infinite-dimensional 
functional spaces. 



II. THE CLASSICAL BRACKETS 

We consider a mechanical system, with coordinates 
(q , ■■■,q n ), described by an action functional S — 
J dt C(q l ,q l ,t), where dot denotes a time derivative. We 
assume that the Lagrangian is a polynomial of second de- 
gree in the velocities q l , and that its Hessian d 2 C/dq l dq> 
is a constant, non-degenerate matrix My. We imagine 
that the system is in contact with a heat bath, and we 
assume that the influence of the heat bath can be de- 
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scribed, effectively, by a mean force, characterized by a 
bounded memory function ^{t — t'), and a random force 
F{ (t) . The time evolution of the system is then described 
by the following equation of Langevin type: 

<)S + f dt'^ ij (t-t / )q 3 (t')=F i (t) . (1) 



Here, 8S/8q l (t) denotes the functional derivative of the 
action S: 
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The original form of the Langevin Equation results from 
the singular limit, where ^(t — t') approaches Jij8(t—t'). 



Mimicking the procedure found by Peierls, to compute 
the Poisson Brackets of a conservative Lagrangian sys- 
tem 0, one can consider the effect, on the system evolu- 
tion, of a small disturbance, produced by an infinitesimal 
change 5S in the form of the action. We consider changes 
of the form 5S = eA, where e is an infinitesimal constant 
and A is a local functional of the trajectory q 3 (t), taken 
from a finite time interval. The small disturbance causes 
an infinitesimal shift 8 A q 3 (t) in the trajectory q 3 \t), and 
it is easy to see that, to first order in e, 8 A q 3 satisfies the 
following linear integro-diffcrcntial equation: 



{L5 A qUt) = - / dt' 
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dt'^t-t') 6 A qi(t')=e— 



(3) 



where it is understood that all functional derivatives are 
evaluated along the undisturbed trajectory. When writ- 
ing the above Equation, we have also assumed that the 
random force does not undergo any variation, to first or- 
der in e. We point out that, by virtue of our assumptions 
on the Lagrangian, the coefficients of Eq. @ depend only 
on the coordinates q l (t) and the velocities q l (t) of the 
undisturbed trajectory, while they are independent of the 
accelerations. This is reassuring, because, by virtue of 
the random external force, the classical trajectories pos- 
sess, in general, smooth velocities, while the acceleration 
does not exist, in the ordinary sense of time-derivatives 
of the velocity @. 
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Since the disturbance A is localized in a finite time 
interval, it makes sense to consider the solution 8 A q 3 (t) 
of Eq.@ satisfying retarded boundary conditions: 



S A q 3 (t) — at early times 



(4) 



The non-degeneracy condition for the Hessian Mij of the 
Lagrangian, ensures that 8~ A ~q 3 (t) exists and is unique. 
We consider also the advanced solution 8\q 3 (t): 



8\q 3 {t) =0 at late times 
of the adjoint equation of Eq. (j3J) : 



(5) 
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(6) 
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where the superscript T stands for transpose (the trans- 
pose coincides with the adjoint, because we are in the real 
field). If B is another functional of the trajectory, with 
support in a finite time interval, we define the bracket 
{A 7 B\ as the following expression, involving the quanti- 
ties Siq 3 (t): 

{A,B} :=-J dt^iS+q'it) S A q\t)) . (7) 



It is immediate to verify that the bracket is bilinear and 



satisfies the Leibniz rule: 

{AB,C} = {A,C}B + A{B,C} , (8) 



{A, BC} = {A, B}C + B{A, C} 



(9) 



To verify that the bracket Eq.JTJ is also antisymmetric 
and that it satisfies the Jacobi identity, it is useful to 
reexpress it in terms of the Green functions G ±K '(£, t';q), 
defined so that: 



8iq\t) 



■ J dt' G ±l3 {t, t'-q) 



8A 



(10) 
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The Green functions G ' l3 (t,t') satisfy the following 
boundary conditions: 

G- ij (t,t';q) = , for t <t' , , (11) 
G +ij \t,t';q) = , for t > t' , (12) 

Mm ^— (t,f;q) = T (M- 1 r '■ (13) 

We define now the commutator function G %J (t, t';q): 

GV(t,t';q) :=G+ i i(t,t' ] q)-G- i i(t,t';q) . (14) 

Note that, by virtue of the boundary conditions satis- 
fied by the retarded and the advanced Green functions, 
G %1 (t,t!) and dtC 3 (t,t') are continuous, in the coinci- 
dence time limit, t — > t' . By using G i0 (t,t'), we can 
rewrite the bracket Eq.Q as: 

^=l dt l dt 'm) &j{tAq) ^W)- (15) 

The antisymmetry of the bracket Eq. Q follows from the 
fact that the commutator function G lJ is antisymmetric, 
as a consequence of the following reciprocity relation, sat- 
isfied by the advanced and retarded Green functions: 

G+ i i(t,t';q) = G-J i (t',t;q). (16) 

Before turning to the proof of Eq. i|16[) . it is useful to in- 
troduce the condensed index notation devised by DeWitt 
0| . With this notation, the trajectory q l (t) is just de- 
noted as q % , with the single latin index i playing the role 
of both the discrete index, and the time variable. Conse- 
quently, repeated condensed indices mean a summation 
on the discrete indices as well as a time integration. For 
example, Eq.Q, with the condensed notation, is written 
as: 

Lij S A q J = (S,ij +n i:j )8 A q J = eA„ , (17) 

where commas denote functional differentiation, and 
Kij 8Aq J is a symbolic notation for the integral linear op- 
erator, depending on the memory function, in Eq.©. To 
prove the reciprocity relation Eq. i|16|) , we point out that 
the Green functions satisfy by definition the following 
Equations: 

G^ k = 81 , (L T )ij G+' k = S k . (18) 

Upon multiplying the second of the above Equations by 
G~ d , we obtain: 

G- a (L T ) i:j G +jk = G- u S k = G- kl . (19) 

However, upon using the first of Eq. l|18fl . we can rewrite 
the l.h.s. of the above Equation as: 

G~ u (L T )ij G +jk = G~ u (L)j l G +jk = 8) G +3k = G +lk . 

(20) 



Upon comparing the r.h.s. of Eq. l|19|l and Ea. (|20|) . the 
reciprocity relation Ea. l|l()|l follows. We can verify now 
the Jacobi identity. Direct evaluation of the quantity 
{{A, B], C} + {{C, A}, B} + {{B, C}, A}, using Eq.© 
shows that: 

{{A, B}, C} + c.p. = A,i B, 3 C, k T^ k , (21) 

where c.p. stands for cyclic permutations of the function- 
als A,B,C. The terms involving second order functional 
derivatives of A, B and C cancel exactly, by virtue of 
the antisymmetry of G %3 . In the above expression, T yfe 
denotes the following quantity, constructed out of func- 
tional derivatives of the commutator function: 

T ijk = QUQjk^ + Qi l G ki ,i +G kl G ij ,i . (22) 

By using the reciprocity relation, the quantity T yfe can 
be written solely in terms of the retarded Green func- 
tion G _u , and its functional derivatives. On the other 
hand, the functional derivatives G~^ k ,i can be computed 
by functionally differentiating the first of Eas. l|18|l : 

l.,,;G ,s ■ l-.j: A -< I). (23) 

Multiplication by G +mi then gives: 

G~ mk ,i = -G +mi L i0 ,i G- ]k = -G-" n L iin i G- jk , (24) 

where in the last passage use has been made again of 
the reciprocity relation. By using this expression into 
Ea. (|22() . it is possible to verify that: 

T i jk = {G -u G - mJG -nk + c . p .) (Lmn;/ _ Lnmil ) i (25) 

where c.p. stands for cyclic permutations of the indices 
ijk. It is easy to check now that T y7c vanishes. Indeed, 
in view of Eq. (|17|) , we see that the quantity between the 
brackets of the r.h.s. is equal to: 

The terms involving third order functional derivatives of 
the action functional cancel each other, because func- 
tional derivatives commute with each other. On the other 
hand, the quantities K%j are independent of the trajec- 
tories q 3 , and hence their functional derivatives vanish 
identically. It follows then that T tjk vanishes, and thus 
the Jacobi identity holds. 

Thus we have shown that it is possible to define a bracket 
on the space of all trajectories. We can evaluate now 
the brackets satisfied by the random force Fi(t). To do 
this, we can use the expression for Fi(t), provided by the 
Langevin Equation, Eq.Q. In this way, we find: 
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{Fi, q k } = {-£,, q\q k } = L i} {qi ,q k } = L l j(G +jk - G~ jk ) = (L - L T ) l3 G+> k + (L T ) y - /.,/; = 



(S,ij S,ji -\-Kij Kji) G — (k k )ij G ^ , 
I 



(27) 



and 



{FuFj} = {-S,i +K ik q k , -S,i +Kji q 1 } = L ik L iX {q\q 1 } = L lk L jt (G +kl - G~ kl ) 



LikLji(G — G fcl ) — Lij Lji — S,ij —S,ji +Kij Kji — (k — n 1 )j 



(28) 
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It is useful to write the above bracket in plain form: 



d/j, 



{F&IFjV)} = ?ZL(t t') + - t) 



dt 



tit 



(29) 



We see from these Equations that, when friction is 
present, the external forces have non-vanishing brackets, 
which implies that they cannot be set to zero. 

Using Ea. (|15|l . it is possible to verify that the equal- 
time brackets for the coordinates q J (t) and the momenta 
Pi(t) = dL/dq t [t) have the familiar canonical form: 



{q i (t),qi(t)} = Q, 



{q i (t),p j (t)} = 6} 



{pi(*),w(*)} = Q 



(30) 



(31) 



(32) 



The verification is similar to the conservative case Q , be- 
cause the memory function contributes to G^(t,t') only 
to order (t — t') 3 . This can be seen by inserting the ex- 
pansions of G ±l i(t, t') in powers of (t — t') into Eas. l(T%)l . 
and exploiting the boundedness of the memory function. 



III. CONCLUSIONS 

We have constructed a set of brackets for a classical 
dissipative system, described by a Langevin Equation, 
with an arbitrary memory function. The brackets sat- 
isfy the usual properties enjoyed by Poisson Brackets of 



Hamiltonian systems. It is worth pointing out the es- 
sential role played, in our treatment, by external random 
forces. When dissipation occurs, they have non-vanishing 
brackets with the system coordinates, and thus cannot be 
consistently set to zero. As a result, our brackets are a 
priori defined on the infinite-dimensional functional space 
of all possible trajectories, accessible to the system under 
the action of arbitrary external forces. However, in the 
absence of friction, when the dynamics is conservative, 
the brackets can be restricted onto the finite-dimensional 
classical phase-space, spanned by the solutions of the 
classical Equations of motion, with no external forces. In 
this case, our construction reproduces Peierls' covariant 
definition of the Poisson Brackets, for dynamical systems 
admitting an action principle. In the context of conser- 
vative systems, the possibility of extending the brackets 
from the phase space to the space of all trajectories, was 
considered some-time ago 0, and our brackets coincide 
with those of Ref.yl, in the absence of friction. 

Quantization can be carried out according to the tra- 
ditional procedure, by replacing the classical brackets 
with commutators. The resulting commutation rules co- 
incide with those that are obtained in standard treat- 
ments of quantum dissipation, by making recourse to a 
microscopic model for the heat bath, after elimination of 
the bath degrees of freedom (see for example Ref. and 
references therein). 
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